On statistics and 1/f noise of molecular random walk in low-density gas 
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The random walk ol test particle in low-density gas is considered basing on approximate coarsened 
version of the collisional representation of the BBGKY equations. The coarsening presumes that 
momentum relaxation rates of the test particle and gas atoms are equal but allows to analyze the 
case when their masses are different. It is shown that both the spectrum exponent and probability 
distribution of 1/f-type diffusivity fluctuations of the test particle essentially depend on ratio of the 
masses, and corresponding distribution of its path is found. 
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I. INTRODUCTION 

In [l|, [4I and Q it was argued (see also (J-Q we argued 
that kinetics of spatially non-uniform gas does not un- 
dergo the Boltzmann equation (or other related and sim- 
ilar equations) even under the low-density or Boltzmann- 
Grad limit ( i^ — > 00 , (5 — > , v5^ = const , with v be- 
ing concentration of gas particles and 5 radius of their 
repulsive interaction). Instead, an infinite chain of ki- 
netic equations appears - for the one-particle distribution 
function (DF) plus infinite chain of specific two-particle, 
three-particle, etc. DFs which represent mean density of 
pair collisions, density of two connected pair collisions, 
and so on. Solution of such chains (as well as their rig- 
orous formulation) will be a great challenge for future 
theory. However, instructive and useful approximate so- 
lutions dan obtained already at present [i|-y, 0, [i3| , first 
of all, for situations when the system (gas or liquid) 
as a whole is almost equilibrium, and the spatial non- 
uniformity is rather of statistical than thermodynamical 
nature. For example, when one considers random walk 
( "Brownian motion" ) of a test (or probe, or marked) par- 
ticle in a fiuid. 

In [l|, 13 and later in [3, [l3| we considered (at differ- 
ent approximations) the case when the test particle, or 
"Brownian particle" (BP), is one of gas atoms. Here we 
consider the case when BP differs from gas atoms (e.g. 
is an impurity atom). 



II. BASIC EQUATIONS 

We star from the Bogolyubov-Born-Green-Kirkwood- 
Yvon (BBGKY) equations [0| for a molecule (molecular- 
size "Brownian particle") dissolved in a fluid [3t3|' 



dt 



= [Hn,Fn]+iy— $fc'(R - r„+i) F„4 

OP Jn+1 






(1) 



where n = 0, 1, . . . , , ^^ ^(r) — V<i>a,b(r) are "atom- 
atom" and "BP-atom" interaction forces, and all desig- 
nations are the same as in [J] or [5| . Since we here again 
are interested in BP's random walk in (quasi-)equilibrium 
fiuid, the initial conditions will be 



F„(< - 0) = A(R - Ro) F,l"')(r(") |R; v) 
xGM(P)n-=iG™(p,) , 



(2) 



where i^n ' (r^^"^ |R; i^) are usual thermodynamically 
equilibrium DF for n atoms in presence of BP occupying 
point R , 

G™(p) = (27rrm)"3/2gxp(_p2/2rm) 

is the Maxwell momentum distribution of a particle with 
mass m , and A(p) is a normalized probability distribu- 
tion like a "smoothed delta- function" . We then want to 
reveal how probability distribution of BP's position R 
evolves after start from A(i? — Hq)!. 



III. COLLISIONAL REPRESENTATION 

If our fiuid is a sufficiently dilute gas, or the 
Boltzmann-Grad gas, then inter-particle interactions can 
be described in terms of well separated collisions. Since 
any "collision" is a finite duration process, while the dis- 
tribution functions (DF) F„ describe instant states of 
the system, we should reformulate DFs and equations 
([T]) in terms of collision processes. 

In principle, this is very simple task. First of all, let us 
consider the (n -I- l)-particle Poisson bracket, [if„,F„], 
describing separate evolution of a group of n fluid atoms 
plus BP, and divide it into two parts, one of which de- 
scribes drift of the center of mass of the group while an- 
other relative displacement of particles inside the group: 
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[iJ„,F„] 



dFr, 9F„ 



9Z„ dQr 



(3) 



We introduced the center of mass positions and velocities, 



Zn 



M-\-nm 



" A/+nm 



and (for n > ) the inner time of the relative motion, 
Qn ■ For example, at h — 1 , 



9Fi dFi ^,, , 



dF\_dF\ 



with pi = ri — R and Mi = Di — V . 

Second, let a given relative disposition of the particles 
can be qualified as a section of collision process (pair col- 
lision at 71 = 1 or sequence of n connected pair col- 
lisions). Then it is reasonable to treat corresponding 
DF Fn as a measure of mean (ensemble average) den- 
sity of such processes, and consider it as a function of 
t, Z„,V„,0„, n relative incoming (m-state) velocities 
(e.g. Uj — Vj — V ) and 3n — 1 input geometric param- 
eters of the process. 

Third, however, such the treatment can not be literal: 
if Fn{t,Qm ■■■) at time moment t represents "density of 
collisions" at this time, then its value at t + At is not 
Fn{t + At,Qn,...) but Fn{t + At,Qn + At,...) , since the 
same collision is represented by different 9„ values at 
different times. 

This statement means that if DF F„ really describe 
collisions as the whole processes then each of EqsH] for 
them (at n > must divide into two equations; 



dt 



-Vn 



dFn 






'Et^/ *a(r,-r„+i)^n+i , (4) 



dFn 
dQn 



, 



(5) 



so that EqlS] must be satisfied "inside collisions" (i.e. 
space-time regions assigned, in relative coordinates, to 
a collision process) while EqjS] "out of collisions" . It is 
natural: since relative movement of colliding particles 
is included, along with their interaction, into collisions, 
it is excluded from description of collision events as the 
wholes (for details and more explanations, see [ll-y, Q)- 
Fourth, the Eqsl5]are the tool for transforming the in- 
teraction integrals in EqlS]into usual Boltzmannian "col- 
lision integrals" : 



dFn _ dFn 



dt 



dZ 



vS. 



n+l ^ n+1 



' ^2-^ ^j n+l-'' n+l 



(6) 



where 3°"^. and Sj^ are the Boltzmann collision opera- 
tors acting on j'-th and A:-th atoms' velocities and on 
BP's and fc-th atom velocities, respectively, and DFs 
F,™j^ describe such {n + l)-particle configurations when 



the "outer" {n + l)-th particle is just on the border of 
the collision. 

Such the border configurations and corresponding DFs 
need in special consideration [i|-y] . If the most (infinite) 
part of the gas stay in thermodynamic equilibrium during 
all (finite-time) evolution of initial (finite-range) distur- 
bance of the equilibrium, then naturally 



Pn+l — Gm{Pn+l) / Fn+1 dpn+1 , 

and Eqsini transform into 

^^" = -Vn • ?§ + ^A^ [Fn+l dpn+1 + 



(7) 



dt 



dZ 



A^ / Fn+lPn+l (8) 






with A"' being the Boltzmann-Lorentz operators. The 
latter are defined by 

S/n+l[Fi-Pr-)GrniPn+l)] = A^" F(...p, ...) 
Sn'+AFi-P-)GmiPn+l)] = A" F(...P...) 

Relations ([7]) and hence EqsJH] do express the "molecu- 
lar chaos" , but with those principal difference from the 
Boltzmann's one that here it concerns velocities of col- 
liding particles only. What is for coordinates of the 
particles, they possess mutual statistical correlations, as 
far strong as strong is non-uniformity of the coordinates 
probability distributions (more comments see PJ, (j, [lyl ) . 
Of course, the transition from EqsH] to EqsjS] should 
be accompanied by corresponding transition from initial 
conditions ([2]) to their collisional representation: 



F„(t = 0,Z,P,p(")) 



O'l 



i^_MR^,nrUr,\ 



M + nm 



i^„(t = 0,fi,r("),P,p("))(ifin 



j=i 



dTj_ 

n 



= 1 






^p-X7E^.-«o X 



3 = 1 



n , n 



i=i 



j=i 



where pj — rj ~ Z and fl is a "collision volume" , i.e. 
volume of a region assigned (in each of the spaces pj = 
Tj — Z) per one collision. Thus under the Boltzmann- 
Grad limit, when a width of the distribution A(p) is krpt 
constant in units of A we can write, obviously. 



i^„(t = 0,Z,P,p(") 
1 



nm\ 



A{Z-Ro) l[GmiPj 



(9) 



This is generalization of initial conditions derived in ^ , 
directly from the Gibbs canonical ensemble, for the case 
when BP is merely one of gas atoms. 

Notice also that under the Boltzmann-Grad limit all 
the variety of the center of mass coordinates Zn can be 
replaced by single common variable Z . 



IV. DIFFUSIVE EVOLUTION AND MORE 
SIMPLIFICATION OF EQUATIONS 

The Eqs|S] still seem very complicated. Therefore we 
would like to further coarsen and simplify them being 
eventually interested in the distributions in configura- 
tional spaces, 



Wnit,Z) 



Fndpi...dpndP 



Since, according to EqsUl 

dt dZ ' 

where Q„ are corresponding probability flows. 



(10) 



Qn{t,Z) 



V„Fndpi...dpndP , 



dQn 

dt 


-A/-/-. 


^I- 


^v,. 


a' + Ta- 



we have to consider also at least evolution equation for 
these flows, 



V„V„F„dpi...dp„(iP+ (11) 
F„+i dpi...dpn+idP 



In order to express right-hand sides here via Q 's and 
Wn 's only and thus close the equations, we need in two 
assumptions. The first is that the conditional velocity 
distributions, Fn/Wn , duffer from the equilibrium ones 
only by non-zero ( .Z -dependent) mean drift velocities, 
Vn = Qn/Wn],^ which is common for all the velocities 
( V and Vj ) under consideration. This hypothesis, in 
turn, can be correct only under the second assumption, 
namely, that all the velocities possess equal relaxation 
rates. Formally, the equality 



Vn 



A" 






= -7K 



approximately takes place, that is P = MV and Pj = 
mvj are left eigenfunctions of operators A'' and A^ , 
respectively, corresponding to the same eigenvalue —7 . 
Then, taking into account also the initial conditions ([9]), 
we can expect that 



/... /'(V„),(V„), Fn dpi...dpndP 



T5, 



M - 



W„ , 



as in equilibrium. Consequently EgsfTT] reduce to 

dQn T dWn 



dt 



AI + nm dZ 



7 Qn+l 



(12) 



and close EgsfTOl 

Of course, on one hand, the second of our assumptions 
looks rather unnatural, since in reality operators A** and 
A^ , may have very different eigenvalues. But, on the 
other hand, the corresponding approximation still allows 
to consider the case of different masses of BP and atoms, 
M ^ m\ and thus obtain a generalization of results of 

i- 

Clearly, our approximation reflects the fact that the 
system "gas -I- BP" is non-equilibrium in statistical (or 
informational) sense only but not in literally thermody- 
namic sense, and therefore its evolution is sooner "dif- 
fusive" than "hydrodynamical" . Nevertheless, if veloci- 
ties of BP and atoms relax (and fluctuate) with differ- 
ent rates then their conditional distribution inside the 
collisional clusters can be significantly non-equilibrium 
(non-Gaussian), and we need in a more complicated ap- 
proximation. 

For the present let us confine ourselves by the EqsfTOl 
and [121 Initial conditions to them, as follow from ([S]), are 

Wn{t = 0, Z) = (1 + ^) A(Z - R,) , (13) 

g„(f = o,z) = , 

where now, of course, we can merely 5{Z) in place of 
A(Z - J2o) . 



V. ANALYSIS OF SHORTENED EQUATIONS 

As usually, let us make the Laplace and Fourier trans- 
forms, writing 

F{p,k) = I dt e-P* jdZ e'^-^ F[t,Z) , 

and introduce designations 

V M V TTl 7 7 

Then EqsHO] and [12] together with EaiT3lvield 

1 _I_ rvji 7 K" 

Wn{p, k) = -^— + - . Q„(p, k) , (14) 



1 + 



X 



1 + an 



P P 



Qn{p,k)^'-^-^Qn+i{p,k) (15) 



p^ p 



According to characteristic structure of the chain of equa- 
tions (jlSp , its solution can be written as sum of infinite 



iteration series. This yields 

^2 oa s 



Qr 



ikVr 



En 



s—n r—n 

CXD S 



7 1 ~\- ar 



1 + an X v-^ T-r 

" s^n r—n 



P'y ^-^ \ V X -\-l ^ ar 

7 1 + ar 



p X + 1 + ar 



(16) 



This is generahzation of of the series from [3] . 
Next, let us introduce quantities 



1 X 1 vlk'^ 

a= hn, , c— 1 Vn = — 5 — I hn 

a a a p'' a 



v^k^ 1 



and rewrite expression (|16p as follows, 



I + an X sr-^ I 7 

^" = + -E hi 



p 



7 V r(a + s) 



7 ^ V P/ r(a)(s-l) 



B{c,n) = 



1 + an _ X^ {-i/pY 



7^r(a)(s-l)!yo 



^a+s-l g-. ^^ ^ 



1 + an 



Jo 

1 / ^ T^-T dt 



a Jo {1 + -it/pY+^ 



Hence, 



pWn 



1 / ^ ^^dt (17) 



1 + an a 7o (1 + ^t/p) i+"+ 

which is direct generalization of formula (27) from |3l]. 

VI. LARGE-SCALE ASYMPTOTIC 

As we already know Jl'-'S], most important statistical 
characteristics of the BP's random walk what distinguish 
it from the usual Ornstein-Uhlenbeck process or similar 
random processes) are its fourth-order path statistical 
moment or cumulant and long-range asymptotic of the 
path distribution. The latter is defined by 

w„(p,fc)= iimc'w^„(e'p,efc), 

g„(p,fc)= lmi^Q„(e2p,efc) 

(for details see Q). We just foreknow that our random 
walk behaves similar;y to diffusion processes (like the 
Wiener process) but with randomly varying diffusivity 
in place of a constant one. Taking into account results of 
[3| , it is reasonable to interpret the asymptotic in terms of 
slow (scaleless) diffusivity fluctuations. Correspondingly, 
let us write 



iy„ = (1 + an) 



Un{D)dD 
p + Dk'^ 



(18) 



where C/„(-D) is effective normalized probability distri- 
bution of diffusivity (at n = it represents BP's diffu- 
sivity while at n > characterizes spreading of densities 
of two- and multi-particle collisional events). 

In the mentioned limit the expression (J17p turns to 



W,. 



1 -|- an 
or equivalently 



1 

-+n- 
a 



1 



exp(-^y) 
pJo (1 + 2/) ^+"+2 

exp {—Dak^r) dr 

(1+Pt)i/"+»+2 



dy 







1 + an 
Combining it with identity 

l/(l + j/)''+i = / 2;''e-(i+«)^da;/r(6- 



(19) 



with b = 1/a + n 



1 , it is easy to obtain 
1 f°° x^ exp (- x) dx 



1 + an r(5) Jo P^ + Dak'^ 
Comparison of this formula and ()18p yields 

1 



Un{D) = 



T{b)D \ D 



exp 



Dc 
D 



-TT (20) 



with h = 1/a + n+l = ii-+n + l. 

Hence, the effective long-range diffusivities are random 
quantities with very bad statistics. In particular, its their 
most probable (m.p.) and mean values essentially differ 
one from another: 



(m. p. D)r 



Dr, 



2 + M/m + n 



(D) 



Dr, 



M/m 



Notice that (!?)„ == Db . 

Making in ([T5| the inverse Laplace and Fourier trans- 
forms and substituting (j20l) . in the space-time represen- 
tation we have 



Wn{t,R) 

00 

{^■RDty^/^ 



lim s'^Wnis^t.sR) 



(21) 



exp 



R^ \ 



ADtJ 

T{b + d/2) 



Un{D)dD 



(AnDaty/'^rib) {l + R^/ADat)''+d/2 ' 

n , R= Z 



where b — 1/a + 1 + n = M/m + 1 
d is the space dimension (d ~ 3). 



and 



VII. STATISTICAL MOMENTS AND 1/F NOISE 

Expansion of the Wn(p,k) into series over k^ gives 
Laplace transforms of equilibrium statistical moments of 
the distributions Wn{t, R) : 



Wnip,k) ^ 1 , ^ (-fc2) 

1 + an 



2\s roo 






t (2^)! Jo 



-pt /n.2s 



(i?^^(t))„dt (22) 



Here and be;ow R{t) means projection of vector R{t) 
onto arbitrary fixwd axis and 

{R^'{t)),,^ J R^'Wn{t,R)dR 

Direct comparison of (f22|) with (fT7|) shows that 

(2.)! 



e-P'{R'%t))ndt 



l(vj 

(s — 1)! p \p^ 



^ (1 - x)i/°+"-i 



In- 



1 



1-x 



X 
s-1 



(23) 



dx 



Considering the long-time asymptotic of these expres- 
sions, i.e. the hmit p/7 — > , we have, firstly. 



-P'{R^it))^dt 



2Vi 



2Dh 



(1 + an)p'^^ (1 + an)p'^ 



(24) 



which corresponds to the usual diffusion-law asymptotic 
{R'^{t))n^2Dbt/{l + an) at -/t ^ 00 . 

Secondly, an asymptotic behavior of the fourth-order 
moment crucially depends on sign of M/m + n — I . 
Namely, if this quantity < 1 , that is n = and 
M/m < 1 , then 



'P'{R%t))odt 



2ADl 

\z~l 



Ci-] ,(25) 



--/"^i^l^.«- 



1 



dx 



m J 

■KZ 



At n = and M/m = 1 , we come to result of [3|, 



(26) 



-P'{R\t))odt 



2ADi 



In^ 

70 P"" P 

And at M/m > 1 under n — we find 

M 



-P'{R\t))odt 



2ADl 

p3 



M -m 



(27) 



(28) 



These three asymptotics should be compared with 
2ADl/p^ which is asymptotic of (i?'*(i))o for the ideal 
Gaussian random walk. 

What is for the higher-order moments at 71 = and 
any moments at n > , they are presented by table 



P'{R^^t))^dt 



(2s)! 1 fvls 

^ > ■-[^] X (29) 



b+l 



(?)' 



Jo 3;''+i 



In 2: 



In^ 



(s — 1)! p yp^ 

s~l 

dx if s > 6 + 1 



(?)• 



r{s)r(b+i-s) 
r(f,+i) 



if s = b + l 
if s <b+l 



with b — 1/a + n . The first row here shows that high 
enough moments are determined by not only the char- 
acteristic diffusion lengths ^/TlT/ji and \/2Dat but also 
by free-flight lengths VqI and vq . 



In the time domain, the asymptotic behavior of 
fourthOorder cumulamt of BP's path what corresponds 
to (ES-lMl) is 



{RHt)),-HR\t))l 



(30) 



^3(2^.^.)^(,.)--/™,f^ 
-^ i{2Dltf ln(7i) 
-^ i{2Dltf "^ 



M — m 



M ^ 
if — < 1 
m 

M ^ 
if — = 1 
m 

M ^ 
if — > 1 
m 



Interpreting these asymptotics as manifestation of low- 
frequency fluctuations of BP's diffusivity, we see that 
in any case (i.e. at any mass ratio M/m ) that are 
non-ergodic fluctuations represented by formally non- 
stationary random processes, in the sense explained in 
1]. Corresponding effective spectral densities of the dif- 
fusivity fluctuations at frequencies <C 7 are 



M/m 



2TTDI 



(31) 



M 

(3-AI)(2-M)eos[f(T3f)] ^' m< 






nDl 



m, 
M ~m 



M ^ 
if — = 1 
m 

M 
if — >1 

m 



The latter expression means "static" fluctuations instead 
of a random process, pe, this is consequence of ap- 
proximate character of our consideration, while the ex- 
act one would change the factor 5{uj) to something like 
w-iln'^(7/w). 

At last, consider asymptotic of sufficiently high-order 
moments. Namely, at s > b = 1/a+n+l = M/m,+n+l 
from (l23t we have 



e-P'{R^'{t))ndt 



(2s)! 1 fPq 

(s- 1)! p\p 



b / 2\ ^^b 



{l-X 



ifc-2 



In- 



X 

s-l 



dx 



and correspondingly 

In 



(32) 



(2s)! 



il~x) 



b-2 



1 



(s - 1)! r(2s - b) J a x'' p l~x 

The coefficient here can be estimated as 



dx 



{Datf (vit^y-b 



< 



(33) 



< 



(2s)! 



l)!r(2s- 



y^-'-'{l + y)'e-^''-'^ydy 



Hence, for fixed b = M /m 
we can write 



1 and large enough t 



IX. RESUME 



lim {R^^{t))]!^ 



< 



vit' 



M/m- 



Tliis means that the distributions Wn{tlR) are sharply 
cut off at characteristic ballistic-flight lengths \R\ ^ 
t^T jijs/l + nm) . In particular, the BP's path distribu- 
tion W^itlR) is cut off at \R\ ~ V^t . 



VIII. VIRIAL RELATIONS 

In conclusion let us concentrate on W^'s dependence 
on the gas density v . Obviously, it is the same as de- 
pendence on 7 , since 7 ex i/ is the only z^-dependent 
parameter of our model, and therefore ^d/ d^ — vd/dv . 

From Eq[T7] it follows that 



Wn+i 



l + a{n + l) 



1 



I + p/7 d 



I/a + n + 1 dv 



Wn 



In the long-range limit it turns into 



W 



n+l 



l4-a(n + I) 



I 



I 



I/a 



T'ai; 



w„ 



I 



an 



(34) 



(35) 



At that, in application to (I^Ul) . we cam make change 
vd/dv = -Dbd/dDb - Dad/dDa . 

In essence, formulas (p4|) - (|35l) represent straight anal- 
ogy of the "virial relations" investigated in 4-8] on rigor- 
ous and most general level. Special treatment of relations 
like (1351) will be done elsewhere. 



In this paper, following [l| and [3|, we continued ap- 
proximate analysis of general equations of the collisional 
approximation to kinetics of spatially non-uniform gas 
lli (see also Sec.III above and [2], [|^ and 0). We ex- 
tended method of [3] ('"diffusive approximation") to the 
case when the test "Brownian" particle (BP) has a mass 
different from mass of the gas atoms, although posse 
the same friction, or momentum relaxation rate, as the 
atoms. By the example of this specific but interesting sit- 
uation we demonstrated that both the spectrum of low- 
frequency 1/f-type fluctuations in BP's diffusivity and 
effective (time-smoothed) probability distribution of the 
diffusivity essentially depend on the mass ratio of BP and 
atoms. This means that probability distribution of BP's 
path depends on this ratio, at that always possessing es- 
sentially non-Gaussian diffusive long-range asymptotic. 

In the case we considered under the approximation we 
used the diffusivity I/f-type, or "flicker", noise has the 
exponent equal to or greater than unit. It remains un- 
clear whether molecular random walk in a gas (or in a 
liquid) can have exponents less than unit or at least "a 
little less" as in the phenomenological theory suggested in 
[l3,ll3| for charge transport (see also [21 )• To answer this 
question and consider more general situations (first of all 
the case of different relaxation rates) we should leave the 
just exploited approximation (since it presumes presence 
of quasi-equilibrium inside the "collisional clusters) and 
start again from the equations ([5]) [ll or even from for- 
mally exact equations of molecular Brownian motion (see 
Sec.II and ^M)- 
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